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Within the a — uo model of coupled nucleon-meson systems, a generalized relativistic Lenard- 
Balescu-equation is presented resulting from a relativistic random phase approximation (RRPA). 
This provides a systematic derivation of relativistic transport equations in the frame of nonequilib- 
rium Green's function technique including medium effects as well as flucuation effects. It contains all 
possible processes due to one meson exchange and special attention is kept to the off-shell character 
of the particles. As a new feature of many particle effects, processes are possible which can be inter- 
preted as particle creation and annihilation due to in-medium one meson exchange. In-medium cross 
sections are obtained from the generalized derivation of collision integrals, which possess complete 
^ , crossing symmetries. 
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> ' I. INTRODUCTION 

\o ' 

Since Walecka [?] has established his model of meson exchange already 20 years ago, field theoretical models 
were successfully applied for the description of heavy-ion collisions at intermediate energies [?,?,?,?]. Thereby the 
derivation of transport equations starting from a microscopic model is paid great attention [?,?,?]. 

The most powerful and appropriate method to describe quantum many particle systems under extreme conditions 
like nuclear matter or condensed stellar objects are the real-time- Green functions technique. Many excellent reviews 
are published concerning this method during the past decade [?,?,?,?,?]. Moreover there are many investigations in 
classical relativistic treatment [?,?]. Several papers are devoted the formulation of transport equations in quantum 
-h field theory [?,?,?,?]. 

In this paper we like to present briefly a derivation of relativistic transport equations for the a — u> model. The 
formal development of relativistic transport models is well known and documented in the literature [?,?]. Here 
• • ■ special attention is kept to the description of fluctuations and the influence of medium effects on the kinetic equation. 
Additional reaction channels are opened by the influence of many particle effects, which would be forbidden in 
uncorrelated systems. Further, we derive a decoupling between nucleon and meson-equations, which result in an 
effective squared one-meson exchange potential containing no mixed terms of meson contributions. This is established 
by the use of generalized optical theorems. 

The outline of the paper is the following. In chapter 2 we review the fundamental equations of the Yukawa 
Lagrangian and introduce the nonequilibrium Green's function technique based on the very fundamental principle of 
weakening of initial correlation. The structure of equations are developed in the Schwinger formalism, generalized to 
the set of four nonequilibrium Green's functions as it was similar done by Bezzcridcs [?] and recently by Davis [?]. 
The spectral information and therefore the quasiparticle properties are discussed in chapter 3 which establish a quite 
natural generalization of Bruckner theory by the aspect of complete particle - antiparticlc symmetry of interacting 
quasiparticles [?]. 

To set up a consistent kinetic theory an equation for the Wigner distribution function should be available, which 
includes all features of relativistic quasiparticle behaviour like Pauli-blocking, screening and relativistic effects, such as 
scattering between particles and antiparticles and pair creation and destruction processes. The consequent treatment 
of many particle effects adopted in section 4 as it was done in the nonrelativistic case [?]. The generalization of 
collision integrals include the effect of density fluctuations by an dynamical meson exchange potential and describes 
completely the particle-antiparticle symmetry. Particularly, it contains Pauli blocking and inelastic processes by the 
medium. A dynamical potential enters the equation due to density fluctuations and the equation can be considered 
therefore as a generalization of quantum mechanical Boltzmann equation [?] to a Lenard-Balescu-type one. This is 
the main result of the present paper and should be the starting point to kinetic description of hot nuclear matter 
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instead of the relativistic Vlassow treatment [?] . Similar equations have been obtained in T - Matrix approximation 
but without pair creation and dynamically screening by De Boer [?], De Groot [?], Botermans [?] and Malfhet [?]. 

From the found collision integrals we derive in-medium cross-sections, which differ from ordinary Born approxima- 
tion by two facts. Firstly, no mixed coupling terms occur between different kinds of meson contributions. This is due 
to the more general approximation of self energies instead of the decoupling normally used. Secondly, the medium 
effects the cross section by an energy shift from the vector-meson exchange and the renormalized effective mass by 
the scalar-meson contribution. 



II. MODEL AND BASIC EQUATIONS 

The Yukawa Lagrangian, which couples scalar and vector meson fields to fermion fields describes the simplest form 
of a model for nuclear matter [?] 

l= -*(-«y*Sp + «o)* + 

+ \d^d^>$> - ^m 2 s <P 2 + : ff s *$* : 

- \f^F^ + ~m 2 v <pxv x + : .g„*7A^ A : ■ (1) 

Here m° s , m° are the bare scalar and vector meson mass and kq is the bare nucleon one respectively. Further, :: denotes 
the normal ordering procedure to ensure the renormalization of baryon density [?] . After renormalization of the masses 
the coupling constants have to be chosen in such a way, that on one hand, the equilibrium ground state properties 
can be fitted, and on the other hand, the scattering data can be reproduced [?,?]. Because of the renormalizability 
of the model Lagrangian we can choose for the occuring divergent vacuum terms a standard procedure resulting in 
physical masses. This will be shown in the Appendix for our used selfconsistent Hartree-Fock approximation. In the 
following we use therefore already the renormalized fields and masses. 
From (Q) the equations of motion read 

(□ + m 2 s ) $ = g s : ^ : = fa (2) 
(□ + m 2 ) ip v = g v : *7„# : = fa v (3) 

(*7 <9 M - k) * = -g s :V$:+g v : 7 A ^ A * : = V (4) 

With the help of the Green's function G„, Gq for the free meson equation (||, ||) the mesonic degree of freedom in 
eq. (0) can be eliminated [?,?,?] 

(i 7 -«)*„= [g 2 (^U (GoT p (ijse - g 2 G s ) (5) 

Here a greek letter indicates the spinor index. Recalling that G S Q — (— k 2 + m 2 s )~ l one may construct from (||) a 4- 
vector potential, which has the structure U s — U v . Therefore in any perturbation expression, where the potential enters 
squared, one gets mixed coupling terms between the vector and scalar mesons. This usually used decoupling overlooks 
the problem of meson initial correlations, which is assumed to be zero when the mesonic degrees are eliminated by 
inverting their differential equation into an integral one, which was done using the free Green's function. Otherwise 
some initial terms have to occur. 

We will proceed in an other way and develop a systematic many particle theory, which yields expressions without 
any mixed coupling terms with respect to the different meson contributions. 

The physical properties of the system can be described by means of Green's functions or correlation functions. For 
nucleon system we define 

iG = (T^W) = iG++ iG< = - = -iG+- 

iG = - (TJSfW) = iG— iG> = (<F*F> = iG-+. (6) 
and for the meson fields we introduce the correlation functions 
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id = (Ttpip') - (ip) {ip') = id ++ 
id= — {T vv i) + (ip) (ip') = id — 
id < = - (f'(p) + (tp') (ip) = 

id£> = {<p^)-{<p){<p')=id- + . (7) 

Here we used the convenient matrix - notation where the ± terms of the 2x2 correlation matrix will be signed by 
latin letters. For the reason of legibility we write down only the scalar meson equations in the following. The vector 
meson Green function carry spin indices. All final results will be presented for both, vector and scalar mesons. 

Using the definition (^) and (0) the equations of motion for the different Green's functions can be derived. The 
equation of motion for the causal one reads e.g. 

(i^d^-^G = 8(1-1') + \{Tj^) (8) 

(n+m s )d = -6(1-1')- -{Tj^ip') (9) 

i 

with jtp and j$ from eq. (j^). Up to now we have not specified the meaning of the average. Solutions, which coincide 
to some special averages, have to be chosen by appropriate boundary conditions. For equilibrium, which corresponds 
to the average with the grand canonical equilibrium density operator, the famous KMS condition [?,?] holds 

G > \ t=o =±e^G < \ tl= _ t0 . (10) 

In the nonequilibrium situation and for real time Green's functions this condition is not valid. Especially G < and G > 
are independent functions. One powerful possibility to derive real time Green functions is the condition of weakening 
of initial correlation [?], which is to be expressed in systems with finite densities. If the condition for the correlation 
time T corr and the mean free collision time t co u 

7~corr ^~coll 

is valid, we have a condition for the right sides of (||) and ([)[) 

lim (j v 4>) = (V^){ 5s ($)-.g,7 A (^A)}. (11) 

This is an asymptotic condition, which breaks the time symmetry and provides irreversible evolution in nonequilibrium 
systems. 

In the next step the correlated self energy E c is introduced formally by subtracting the mean field parts [Til) from 
the right side of eq. (||) 

- (Tj^) - lim (Tj v ip) = f E C G. (12) 



The way of integration c has to be determined in such a way that ( |ll| ) is fulfilled. This can be found by 

+ OG 

f dlZ(l,l)G(l,l') = J dl{E(l,l)G(l,l')-S < (l,l)G > (l,l')} 

c — OO 

= £ / S +d G d+ . (13) 



It is easy to see that the boundary condition is fulfilled, since the contribution (|l3| ) vanishes in the limit t' x = t^ 1 — > — oo. 
For the case t\ < t[ (and vice versa) we can write e.g. 
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dl{S(l,T)G(i,l')-S < (l,l)G > (l,l')} 



E ; (l. IKV ( I . I ' ) - r /'e il. DG + 

— oo 
oo 



+ / E<(1,1)G > (1,1')- / S<(1,1)G > (1,1 / ) 

for ti->ti->t = -oo (14) 



If we split the last integral on the right of (13) into two parts according to 

dti 



p + OD 






/ dh = 






/ — OO 


/ — OO 





a contour of time integration follows which is equal to the Keldysh-contour [?,?]. 

Now we are able to enclose the equation of motion (||) resulting in the nonequilibrium Dyson equation as matrix 
equation (^]) in the following form 

oo 

(iys,, - «) g 6c (ii') = <5 6c (n') + /" d T (i T) Ga C (T i') (15) 

—00 

and for the meson correlation functions one has 

(□ + to 2 ) 4c(H') - S bc (l - 1') + J d2U bd (12)d dc (21 r ). (16) 

d 

Equations (eqll) and ( |l6| ) are matrix equations, where latin letters remark (+,-) in the following. To study the 
structure of the self energy S and the polarization function II explicitly we want to use a technique developed by 
Schwinger [?] . Therefore we introduce an infinitesimal meson generating flux jq, for scalar mesons and j\ for vector 
mesons, as a new type of interaction 

hnt = -]^-h^ X - (17) 

In our presented formalism it is not necessary to introduce a generating functional for nuclcons. It turns out that the 
infinitesimal interaction (|l7|) is sufficient to obtain the Kadanoff-Baym [?] equations. Therefore they are valid for any 
density or correlations in the system. 

At this point it has to be remarked, that we suppress the nondiagonal terms of the vector mesons by choosing the 
generating functional in diagonal form. This is possible, because they must not contribute to physical observables 
such as S-matrices as a result of coupling to the baryon conserving flux. The influence to non - observable quantities 
such as self energy should be considered in principle. But they are neglected here for the reason of simplicity, which 
is in agreement with the treatment of Horowitz and Serot [?] . 

Now we introduce an interaction picture in respect to the infinitesimal interaction (|i"7[). Further, we distinguish 
between the upper and lower branch [?] of the contour because in such a way we can find relations for the time 
evolution operator of this (infinitesimal) interaction on the Keldysh contour, which reads 



S , c = T c exp|-ic J (-j*$-^^) c | 



(18) 



Thus special relations can be established by variational technique 

^r- — -ibT b ip' S b S bc , (19) 

where j stand for j$ or j\. Latin letters b, c remark (+,-) respectivly. In a straight forward manner one can express all 
correlation functions ([?]) through variations with respect to j [?]. The causal Green's function e.g. can be expressed 
as 
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^y ( nn >+ = (T+fO-M^')) ~ < >+< v(i') >+= 
<u+(i') 



(20) 



With some manipulations the introduced polarization function (16) of mesons can be derived from (^|) in the form 

n M (12) = -g 2 &]T J Tr[G ba (lT)r a7M (1T2) Ga b (T 1+)} (21) 

aa 

where we have introduced the vertex function 



Taad 112 = - 



1 SG~l a (l l) 

5 <5 < 93(2) 



(l-l) ^(1-2) + ! 



'£gg(l l) 

3 <5 < v?(2) > d ' 



(22) 



The second line of (|2^) is quite easy to verify by means of (|l5|). Concerning the nucleons we find that the right hand 
side of equation (pi) could be expressed by the help of variations in such a way that the equation rtq) takes the form 



[irf8p -n Q +g s <cj)> b -g v ((f X ) 7 A ] G hc (ll') = 5 6c (ll') 
>9s S bg v 6 \ 



(23) 

After the introduction of the same vertex function (p2|), the structure of the self energy introduced in ( |l5| ) can be 
written finally as 



s ba (i 1) = -j2 



bg 2 G ba {ll)T a - ad {l 1 2)^(21) 



(24) 



Equations ([L5|) , (|iq ) , (g_TJ) , (|22| ) and ( P4[ ) form a complete quantum statistical description of the many particle system 
in nonequilibrium and serves as a starting point for any further approximated treatment. We have to point out, that 
the different contributions of mesons are linearly additive in the self energies and therefore quadratic additive in the 
coupling constants as it is to be seen from (p3|). This differs from other treatments, where on the stage of eq. (2) the 
different parts are decoupled approximately yielding a linear additive coupling. Since we have not used any restriction 
or approximation in deriving the general structure of (23), it has to be considered as a more general treatment. The 
physical reason is, that during the derivation of the normal simple additive coupling initial correlations of different 
meson are neglected. By the sceme proposed here this is done carefully resulting in a quadratic additive behavior. 



III. SPECTRAL INFORMATION 

Because a complete description of a nonequilibrium situation demands two independent correlation functions instead 
of one, as it is the case in equilibrium, we have some relations between the four Green functions 



G r -G > =G a -G< =G 
G r + G< =G a + G> = G, 

where we introduce the retarded and advanced functions [?] 

G r ' a = G 5 ± 9 (± (x - x' )) {G> -G<}. 

This function determines the spectral information [?] and we get the equation for them from ( |l5| ) 

K)G r = <5(1 1')+ / dTs r (lT)G r (Tl'). 



(25) 
(26) 

(27) 

After the Wigner transformation, the variables are splitted in macroscopic and microscopic ones and it is assumed [?] 

d d 



£ r (*,p)» 



Z r (X,p) 



> 



(28) 
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This means that Ax M Ap^ 3> 1 where the characteristic length Ap, at which the self energy varies in four-momentum 
space, corresponds to the inverse space-time interaction range. Therefore the approximation (|2^) demands the 
shortness of the space-time interaction range when compared with the systems space-time inhomogeneity scale [?] . 

All this is the physical background, when we applied the so called gradient expansion at this place, which yields 
now for ( p7| ) 

{<fp ll -K-Yr(pX))Q r ( P X) = l. (29) 

To make further progress, we write the self energy in characteristic parts in terms of the 7— matrices, which can be 
verified by general considerations [?,?] 

E r = 7 M S r + IYj r (30) 

Further it is useful to split up the self energy in real i?eS and imaginary T parts, between which exist the dispersion 
relation 

MT(M = pf^M. (31) 
J 2n uj — uj 

Introducing the following medium dressed variables 

M = k + ReHj 

r 



G= (70*06 + -A (7P + M) (32) 
one obtains the complete spectral function from (129]) as 



2G 

A = i{G r (uj + is)-G a (uj~ie)}= ( 7^P M + IM ) g . (33) 

/ 'P 2 — M 2 ! + G 2 



This derivation underlines the correct expression using £7° instead of e in the denominator as it was pointed out by 
P. Henning [?] . In the case of vanishing damping T we get the spectral function, which determines the quasiparticle 
energies in the quasiparticle picture 

A = 2tt {y*P M + / M} S (p 2 - M 2 ) sgn (p ) . (34) 

The 5 Function in ([mJ) represents the defining equation for the quasi-particle-energies 



E 1/2 = ReZ r ± V (p/ - ReUjY + {k + Re^Y \ u=Ei/a (35) 

This is a nonlinear equation by the many particle influence and shows immediately the particle/ antiparticle symmetry. 
Therefore the single particle and antiparticle states become dressed by the medium and yield a natural generalization 
of the Dirac-Bruckner Theory. It can be seen that ( |35|) is determined only by the retarded self energy, which yields 
a mass off-shell behavior of the quasiparticles. This is important to note for the discussion in chapter 4. 

Let us now consider the meson equations, which read for the retarded Green's function written in operator notation, 
where integration about inner variables is assumed 

(□ +m 2 ) (-d r ) = 1 - U r d r . (36) 

Moreover we have for the correlation function 

(□ + m 2 )(-d)< =-U r d< - U<d a . (37) 



Combining (36) and (|37j) we obtain the following important relation 

-d< =d r n< d a (38) 
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which presents the generalized optical theorem [?,?]. Following the same arguments as deriving ( j29| ) one gets for the 
product d r d a in gradient expansion 



1 



d r d a = 

Here the quasi particle energy of mesons is introduced by 

E\ = p 2 + m 2 - Rell. 
The spectral function of mesons in the case of nearly vanishing damping reads 

B = i(d r — d a ) — ► 2ttS (uj 2 - E 2 d ) sgn(w) for T -> 0. 



(39) 



(40) 



(41) 



In static case with vanishing damping one obtaines from ( |39| ) the Fourier-transformed Yukawa-potential for the 
meson exchange 



d r d a 



V(r) 



1 



5 4 (4^) 4 



V 2 (p) 



r~ 2 = m 2 - Ren. 



(42) 



When II is determined by ( pl| ) we have the possibility to construct a self consistent system including fluctuation 
phenomena. 



IV. KINETIC EQUATIONS 



The starting point are the generalized KB equations [?] which are exact and read in the space-time notation 



ReG r , G 



, ReG 



(43) 



Here we used operator notations, where the integration is assumed over inner variables and the brackets sign the 
commutator [] and anticommutator {} respectively. 

Using the gradient expansion here we have a different physical content than the one yielding ( p9| ) and (|36|). If we 
suppose here 



G(X,p) » 



d d 



dX» d Ptl 



G(X,p) 



> 



(44) 



we have to demand that AXAp 3> 1. This would be justified in the case of one particle systems only by classical 
description. Because we have a many particle system G contains information averaged over space-time cells, which 
can be in principle smaller than the one determined by the single particle de Broglie wave length [?]. Because in 
nuclear physics time gradients occur, which are in principle nonvanishing, this assumption is the most restrictive one 
and is now under investigation and will soon be published. We have from (^3|) in pujRT variables 



ReG r ~ , G < 



= G< £> - G > E < 



where the bracket means the Poisson-bracket: 



[A > B] = iL a -^t b - L B ik A + ^ rA ^ a v * b - 



(45) 



(46) 



As a first step of approximation we neglect higher vertex corrections in agreement with the gradient expansion and 
derive from ( |22|) for the vertex function 



(47) 



Taadl 112 «5 oS (l-l)<5 od (l-2). 
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In the framework of this relativistic Random Phase Approximation (RRPA) we get from (|l5|) , ( |l6| ) , ( pl| ) , ( p4| ) , ( |38| ) 
and (|47| ) the complete set of equations 



G r ,G < 



(puRT) = G<£> - G > S > 
£<(xA) - -z.g 2 G< (l f) d< (T 1 



d<(ptuRT) = —L^V 2 {p^)U< 

> r > > 

n < (xX)=ig 2 Tr\G < G < 



(48) 



It has to be stressed that the dynamical potential introduced by the optical theorem ( p9| ) now contains an infinitesi- 
mal sum of nucleon fluctuations by the polarization function resulting in an effective meson exchange mass. Therefore 
this approximation can be considered as a modified first Born approximation including collective effects, especially 
density fluctuations. 

> 

The set of equations (|48|a,-d) forms a closed set and determines the correlation function G < and therefore the kinetic 

> 

description proposed we find a connection between G < . 

Therefore we proceed to the kinetic level of description by introducing the generalized distribution 



G > {puRT) = -iA(puRT)( 1 - f(puRT) ) 
G< [puRT) =iA{puRT) f{puRT) 



(49) 



which is until now an exact variable change without approximations. Especially it fulfills the conditions (p3|). 

As far as the generalized distribution should describe physical particles and antiparticles we use the Dirac interpre- 
tation in the quasiparticle picture. There an empty state of particle with negative energy equals to an antiparticle 
state with positive energy and we can write 



f(-u>) = 1 - F(u), 



(50) 



where F signs the antiparticle distribution. In equilibrium this is identical with the fact, that the chemical potential 
has to be chosen with opposite signs for particles and antiparticles, which follows immediately from the conserved 
baryon density [?] . 

With the help of (§2) and (§§) we obtain for G< 



G K 



{ (7°£'i - 7 X Pi + M) f (£i) 5(u- E x ) 
-j°E 2 - 7 : Pi + M) f (-Ea) 8{u + E 2 )}. 



(51) 



The choice of (pfl) coincides with Bezzerides and DeBois [?] if we set F < (lu) — f(u>) and F > (uj) = 1 — F(u>). 



we have to perform the Poisson brackets on the 
re renormalization denominator, which arise from the 



If we now introduce (|50j) and ( pl[ ) into the kinetic equation (4 
left side carefully. After partial integration one can show, that t. 

spectral function (|34j) cancels exactly with the factors following from the gradient expansion on the left side so that 
the drift term takes the form 



d 



we want to restrict to virtual meson exchange and neglect all density terms in the meson Green's function. After 
somewhat extensive but straight forward calculation using the set of equations (Ji^a-d) we finally arrive at the kinetic 
equation for the particle distribution, if we integrate both sides over positive frequencies. The equation for antiparticles 
are obvious by replacing f p < — ► F p . 

We use for shortness the following abbreviation of the dynamically potential ((3^) without static approximation 



V 2 {p,p,pi,pl) = 

(2(pp 1 )(pp 1 ) + 2(pp 1 )(pp 1 ) 

- (pp + K 2 )(piPl + K 2 )V 2 (p - 



- 2 K 2 [(p 1 p 1 ) + (pp)} + 4 K 4 ) V 2 (p - p) 
P) 



(52) 
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As already mentioned following ( |24| ) no mixed terms between g r , g s occur in the cross sections due to the additive 
behaviour of the different self energies. The kinetic equation reads now 

jjjif + VrE Vpf- V P E S7Rf = Tr J [g<Z> - 5 >E<] ^ + 

with the collision integral up to first two order gradient expansions. The second part of the right side of (^) is 
obtained by the second order gradient expansion or equivalently by the expansion of the time retardation up to first 
order. It ensures the complete energy conservation [?], whereas the first part of (|53| ) ensures only the quasiparticle 
energy. Because we restrict here to the case, where the quasiparticle picture is valid, only the first part will be 
discussed. In a later paper we will present the results obtained for the second part. The first part of the right side 
from (|5^) containes 8 processes. The first one, which describes elastic particle particle scattering reads 

d P l dpi (27r)V 4 )( - i + -_ pi _ p) . 1 



(2tt) 3 (2tt) 3 (2tt) 3 y ' ^ r ^ L ' \E X E X EE\ 
V 2 (p,p,p 1 ,p 1 ) {f P f Pl N (f p f Pl ) - f p f Pl N (f p f Pl )} . (54) 

Here and in the following we used the Pauli blocking factors 

N(F) = 1 - F 
N(Ff) = (l-F)(l-f) 
N(FfG) = (l-F)(l-f)(l-G). 

The next two processes are the the crossing symmetric processes in the t and u channel. They are obvious by the 
crossing symmetric relations and describe the elastic particle- antiparticle scattering. The next 5 processes are only 
possible due to the off-shell character of the quasi-particle energies ( |35"| ) 

dp 3 dpi dpi 4 1 
(2tt) 



(2tt) 3 (2tt) 3 (2tt) 3 v ' \E 1 E 1 EE\ 
[ - V 2 (p,p,-pl,pl)S^ (px+Pi +p-p){f P N(f p F p J Pl ) - f p f Pl F Pl } 

- V 2 {p,p,p 1 ,~p 1 )5 {4 ) (~pi - pi + p ~ p) {f p f pi F pi - N (f p f Pl F Pl ) fp} 

- V 2 { Pl p uPu -p)5^) (-p-px +pi-p) {f P f pi F p - N (f p f Pl F p ) f p } 

- V 2 (p, -p, - Pl , -pi)5 (4 ) (-Pi +JH-P-P) {f P F p F pl - F Pl N (f P F p F pl )} 
+ V 2 (p, -p,pi,-pi)<5 (4 ) {-pi -pi-p-p)* 

{fpF P f Pl F Pl - (1 - /„) (1 - F p ) (1 - f Pl ) (1 - Fp,)}] . 



(55) 



Analyzing the energy momentum conserving ^-function on finds in connection with ( p5|) that these processes occur if 
the cm. energy y/s fulfills the condition 

= -2ReTT . (56) 

This means, that in a nucleus system, where the zero part of the selfenergy become negative, new reaction channels 
will be opened by many particle influence. One may think on pseudoscalar and pseudovector meson coupling, where 
the absence of mean field terms of the vector mesons ensures the condition (j56"l). 

Further, it can be seen explicitly that the relativistic treatment yields some forefactors ( |52| ) to the dynamical 
potentials (^9|). The dynamical potentials have to be determined selfconsistently with the kinetic equation by means 
of ( |5li| ) and the polarization function (|48|d). 

Therefore we arrived at a generalized relativistic Lenard-Balescu equation. The dynamical potentials reflect the 
influence of virtual mesons and are written for each process with the influence of the dynamical density fluctuations 
by ©• 

The first 3 terms are presenting the elastic scattering between particles and antiparticles respectively. The next 5 
processes are only possible, if the particles are off shell, which is ensured by the quasiparticle energies shown in (|35|). 







This can be understood as particle creation and destruction by means of virtual meson exchange. It is important to 
recognize, that these processes are special effects of many particle treatment and are forbidden by momentum and 
energy conserving <5-functions in the case of vanishing many particle influence. 

The equation ( p3|) possesses all properties of a selfconsistent quantum-mechanical kinetic equation which describes 
the relativistic behaviour of a meson-nucleon system from the many particle point of view. Particularly, it contains 
Pauli blocking and inelastic processes by the medium. This is the main result of the present paper and should be the 
starting point to kinetic description of hot nuclear matter. 

To complete we give the in medium cross sections obtained for the elastic process (^4|). For shortness we denote 
ReY, r = A and k = k + i?eS[. The proton-proton cross section reads 



|)» + (A(VJ+A)-S 2 + f) 2 



Here s,t,u are the Mandelstam variables. This means that the in medium cross sections are effected by the dressed 
Baryon mass and an additional term proportional to the vector part of the self energy. 

The in-medium proton-antiproton clastic cross sections as well as the discussed medium caused inelastic ones can be 
derived from the corresponding collision integrals ( |54|) and ( |55| ) in the same way. It turns out that they are identical 
to the cross sections one obtaines from ( |5^ ) by using crossing symmetric relations in the t or u channel or to inelastic 
crossing. This means that we end up with an expression for the medium dependent cross section for nucleons, which 
possess crossing symmetries and includs an infinitesimal series of meson interactions by the RPA polarization function 
resulting in an effective meson mass (fIo"|), 

V. CONCLUSION 

We have analysed the structure of equations for a nonequilibrium situation consisting in a relativistic nucleon- 
meson system. New features arise due to the consequent nonequilibrium treatment as pair creation and destruction 
and dynamically interacting meson masses. With one formalism of real time Green function one can determine the 
self energy, the quasiparticle behaviour and kinetic equations of the system and therefore the scattering measures. 
The kinetic equations are derived of Lenard Balescu type by V s approximation which coincide with the RPA screened 
potentials in nonrelativistic treatment. These obtained potentials may serve as a starting point for constructing optical 
potentials [?]. On the other hand in-medium cross sections can be derived from the kinetic equation, which show no 
mixed terms of mesons contributions. This follows immediately from the more general decoupling of equations by the 
self energies. In a forthcoming work we will present the numerical results for the in-medium cross sections derived 
from (|5^) for the now opened inelastic scattering channel. 
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